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Abstract 

We establish a link between the two functional approaches: a mesoscopic field theory developed 
recently by A.Ciach and G.Stell [A. Ciach and G. Stell, J. Mol. Liq. 87 (2000) 253] for the study 
of ionic models and an exact statistical field theory based on the method of collective variables. 



Two rigorous scalar field theories were recently developed to describe the phase equilibria 
in ionic fluids: m the KSSHE (Kac-Siegert-Stratonovich-Hubbard-Edwards) theory and 



the approach BJ3 [which is based on the collective variable (CV) method |4|, |5|. As was 
shown recently 13, y] both of theories are in close relation. Our goal here is to establish a 

nn 

link between the CV approach and the mesoscopic field theory developed in p, |8( for the 
restricted primitive model with additional short-range interactions presented (RPM+SR) or, 
more specifically, to demonstrate how after some approximations in the exact microscopic 
CV action TC[un, p, Q, u, 7] we can arrive at the functional of the grand potential AQ MF [i], 0] 
considered in p, |S| . 

Let us consider a general case of a classical two-component system consisting of TV par- 
ticles among which there exist Ni particles of species 1 and N 2 particles of species 2. The 
pair interaction potential is assumed to be of the following form: 

U a p(r) = <?(r) + v%(t) + v™(r), (1) 

where v^r (r) is the interaction potential between the two additive hard spheres of diameters 
a aa and a pp. We call the two-component hard sphere system a reference system (RS). 
Thermodynamic and structural properties of RS are assumed to be known. v^Jr) is the 
Coulomb potential: v^Jr) = q a qpv c (r), where v c (r) = l/(Dr), D is the dielectric constant, 
hereafter we put D = 1. The solution is made of both positive and negative ions so that the 
electroneutrality is satisfied, Yla=i QaCa = 0, and c a is the concentration of the species a, 
c a = N a /N. The ions of the species a = 1 are characterized by their hard sphere diameter 
(Tii and their electrostatic charge +q and those of species a = 2, characterized by diameter 
(X22, bear opposite charge —zq (q is elementary charge and z is the parameter of charge 
asymmetry), ff^(r) is the potential of the short-range interaction: v^§(r) — v^Jr)+v^Jr), 
where v^Jr) is used to mimic the soft core asymmetric repulsive interaction, v^Jr) is 
assumed to have a Fourier transform; v^Jr) describes a van der Waals-like attraction. 

We consider the grand partition function (GPF) of the system which can be written as 
follows: 
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Here the following notations are used: v a is the dimensionless chemical potential, v a = 
fln a — 3 In A, /i a is the chemical potential of the ath species, (3 is the reciprocal temperature, 



A -1 = {2nm a f3~ l /h 2 ) 1 / 2 is the inverse de Broglie thermal wavelength; (dr) is the element 
of configurational space of the particles. 

Let us introduce operators p k and Qk : P k = S Q Pk,a an d Qv. = ^2 a QaPu.,a, which are 
combinations of the Fourier transforms of the microscopic number density of the species a: 
Pk,o = Yli ex P( — ikr"). In this case a part of the Boltzmann factor in (J2J) which does not 
include the RS interaction can be presented as follows: 
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and we use the notations $>*p-(k) = ^v^--(k) with v^\k) being a Fourier transform of the 
corresponding interaction potential. 

In order to introduce the collective variables (CVs) we use the identity 
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where 5jr[£ — £] denotes the functional delta function. 
Taking into account we can rewrite (J2J) 



(dp) (dQ) (dw) (d7) exp (-W[i/ a , p, Q, w, 7] ) 



(6) 



where the CV action 7i is as follows: 

H[u a ,p,Q,UJ,j] = - y^[$ArAf(fc)pkP-k + $QQ(fe)QkQ-k + 2$AfQ(fc)pkQ-k] 



-i^(o; k p k + 7kQk) - lnE H s[z/ a ; -iw, -ig«7]. 



(7) 
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In (0) CVs pk and Q k describe fluctuations of the total number density and charge density, 
respectively. E#s[z/ a ; — iu, — ig a 7] is the GPF of a two- component system of hard spheres 
with the renormalized local chemical potential 

+ \ E + I E ^ - + (8) 
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The MF approximation of functional (JUJ) is defined by 

z-mfWoI = exp(-H[v a ,p, Q,U),j]), (9) 
where p, Q, lj and 7 are the solutions of the saddle point equations: 

p = (N[u a ; -iQ, -iq a j]) H s, = 0, 

q = -ipv NN (0), 7 = -ipv NQ (0). (10) 

Now we present CVs pk and Qk (^k and 7k) as 

Pk = p5 k + 5pk, <5k = QS k + 5Q k , 
^k = ^<$k + 5^k, 7k = 7#k + ^7k, 

where the quantaties with a bar are given by (|1U|) and S k is the Kronecker symbol. Then we 
write mS# g[j/ a ; — icu, —1507] in the form of the cumulant expansion 

ln5* ff [...] = EH^T^ E M ( > n) (k 1 ,...,k n )6 lkl ...6 lkin 

n>0 i n >0 ki,...,k n 

X ^k ln+1 • • • Suj k J kl+ ... +kn , (11) 

where DJln n \ki, . . . , fc n ) is the nth cumulant defined by 

mt\h, ...,&„) = 5 ln ^r [ "' ] — a7 — l^=o,^ k =o. (12) 

d8^ kl ...d8^ kin d5u kin+i ...d5uj kn 
The expressions for the cumulants (for n < 4) are given in [3|. Substituting (|ll|) - (fT2*|) in 
©-{ZD we obtain 

s a ] J (d5p){d5Q)(d5uj){d5j) exp | - i E [$Jw(*Ofykfy-k 
+2$ J v Q (fc)5p k 5Q_ k + ®QQ(k)5Q k SQ- k + i E^kfyk + ^7k5Qk) 
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Let us make the following approximations. We neglect the fc-dependence of the cumulants 
putting Wln n \ki, . . . , k n ) ~ Wln"\o, ... ,0) and replace the full chemical potentials u a by 
their MF values z/*. We also limit our consideration to the restricted primitive model (z = 1) 
supplemented by the same short-range interactions for both species (RPM+SR) (vf R (r) = 
v 22 R ( r ) = v i2 R ( r ) — v SR (r)). The latter means that $ N N(k) = $vq(^) = and 

$ QQ (k) = <S> c (k). In this case we have for the cumulants [2j 

*m® = d n , = g 2 G n _ 1; mtf = o, 

m<P = g 4 (3G n _ 2 -2G n _ 3 ), (14) 

where G n = G n (0, . . . , 0) is the Fourier transform of the n-particle truncated (connected) 
correlation function [9j] of a one-component hard sphere system with the density p defined 
by (fTUj) . Because both DJl^ and DJl^ are positive and smooth functions in the region 
under consideration [2], we can integrate in (|T3*j) over 5ujy_ and 5^ with the Gaussian density 
measure as basic one. The integration is performed using the Euler equations. We determine 
Su^ and #7^ which provide a maximum for the functional in the exponent of (j!3)l . This leads 
to the expression for 3 

3 = E MF C J (d5p)(d5Q) exp (-H(6p, 5Q)) , (15) 

where 

H(5p, 5Q) = \ J2[a { 2\k)5 Pk 5p^ + af (k)6Q k 6Q- k ] + E ^[ E 

k n>3 «n>0 
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and the following notations are introduced 

a ( p(k) = ^ SR (k) + l/m^\ af(k) = ^ c (k) + l/mf\ (17) 
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Taking into account ()14)1 we rewrite (jlfij) - (|2T)|) as follows: 



(20) 



S = ZmfC / (d<5p)(d5Q*)exp (-H(5p,5Q 



(21) 



where 



H(6p,6Q*) = ±Yl^\k)8p k 8p„ k + a?\k)5Qt8Q\]+J2^^ 

k n>3 
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r(4) 



(22) 



with 5QI = 8Q k /q and S n = G n /(N). 

After taking into account y = (27r)~ 3 f dk and 5ki+...+k n = y f drexpir(k! + . . .+k n ) 
in (|22J1 we arrive at the same expression for the action as that obtained in the mesoscopic 
field theory (see e.g., Eqs. (44)-(48) in |g|), AQ MF [fj,$\ with 7y(k) = 5p k and J(k) = 
In order to demonstrate that the corresponding coefficients of both actions coincide i.e., 



a { 2 \k) = C^{k)=(3v SR {k)+ lo , 2 



(23) 



a?\k) = C^(k)=^(k)+ l2fi , 



(24) 



= 70,3, 4^ = 3 72,1> = 70,4, 4^ = 6 72,2, 4 = ^ Afl 



(25) 



(C$(k), C${k) and l2m,n are the notations used in |8j), we consider coefficients 7 2m ,„ in 
detail. To this end we recall that 72 m ,n denotes the appropriate derivative of the Helmholtz 
free energy of the hard sphere system f h (see e.g.jjj) 

d 2m+n f h 

iQp* 

For example, when the Carnahan-Starling (CS) form of fh is adobted in the local-density 
approximation 

Rf( * a\ P*±±\ f P* + <t> \ , P*-$ , ( p*~A *^ * s (4-3s) 

PfhKP , <t>) = n l0 S n + n l0 S n ~ P + P 



l2m,n - P fiJOmfln*™ \<f>=0,P*=P* 



2 °V 2 / 2 °\ 2 / ' r (1-8) 

we obtain the following explicit expressions for 72 m , n 

l + 4s + 4s 2 -4s 3 + s 4 1 1 

-"••' = (T^w- = ^5- *» = f (26) 

l-5s-20s 2 -4s 3 + 5s 4 -s 5 S 3 1 
7o,3 — T^rr — l2,\ — - — , I 2 ' ) 

2(l-6s + 15s 2 + 52s 3 + 3s 4 -6s 5 + s 6 ) 1 /- Jf\ 

P (1 - s) b p* 3 Sl \ S 2 J 
2 2 

72,2 = — , 74,4 = — • (29) 

p* p* 

In the above equations we use the notations p* = p* + 5p* and s = up* /6 (p* = pa 3 and 
p* are the dimensionless total number density and its MF value respectively). As is seen, 
expressions (|2fi j) - (j29j) confirm the relations given by (j23jl -(|25 jl . 
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